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ABSTRACT 

During flares and coronal mass ejections, energetic electrons from coronal sources typically 
have very long lifetimes compared to the transit times across the systems, suggesting confinement 
in the source region. Particle-in-cell simulations are carried out to explore the mechanisms of 
energetic electron transport from the corona to the chromosphere and possible confinement. We 
set up an initial system of pre-accelerated hot electrons in contact with ambient cold electrons 
along the local magnetic field, and let it evolve over time. Suppression of transport by a nonlinear, 
highly localized electrostatic electric field (in the form of a double layer) is observed after a short 
phase of free-streaming by hot electrons. The double layer (DL) emerges at the contact of the two 
electron populations. It is driven by an ion-electron streaming instability due to the drift of the 
back-streaming return current electrons interacting with the ions. The DL grows over time and 
supports a significant drop in temperature and hence reduces heat flux between the two regions 
that is sustained for the duration of the simulation. This study shows transport suppression 
begins when the energetic electrons start to propagate away from a coronal acceleration site. It 
also implies confinement of energetic electrons with kinetic energies less than the electrostatic 
energy of the DL for the DL lifetime, which is much longer than the electron transit time through 
the source region. 

Subject headings: Sun: corona — Sun: flares — Sun: particle emission 



INTRODUCTION 



In flares and coronal mass ejections, coronal 



electrons can be accelerated above 100 keV (Lin 



et al. 2003 Krucker et al. 2010 2007), i.e., more 



than three orders of magnitudes higher than the 
ambient coronal temperature of 100 eV. The en- 
ergy is ultimately believed to come from mag- 
netic fields through magnetic reconnection, which 
occurs throughout the corona. X-rays and mi- 
crowaves produced by these energetic electrons 
are observed in both the chromosphere and the 
corona. How energetic electrons are transported 
from a coronal acceleration site to the chromo- 
sphere is a key issue in understanding the dynam- 
ics of flares. Recent time-of- flight measurements of 



rival at the chromosphere by electrons of differ- 
ent energies propagating from the coronal acceler- 
ation site. The energy dependence of these time 
delays provides evidence that electrons can prop- 
agate along magnetic field lines by free-streaming. 
However, hard X-ray observations from above- 



the-looptop (ALT) sources (Masuda et al. 1994 
Krucker et al.|[2"010[ 20071 suggest confinement 



hard X-ray emission during solar flares ( Aschwan- 



den et al.nl996|) report relative time delays of ar- 



of the energetic electrons above the bright soft 
X-ray loops for periods several orders of magni- 
tude longer than the (free-streaming) transit time 
across the source. Increased detection sensitivity 
indicates ALT hard X-ray emission is likely a com- 
mon feature of all flares ( Petrosian et al.|2002 ) al- 
though there have been fewer detections of ALT 
emission than chromospheric emission because the 
higher density, and hence emission intensity, in 
the chromosphere overshadows ALT emission. In 



this case, there must be mechanisms that inhibit 
the transport of electron heat flux. It is puzzling 
why some electrons freely propagate to the chro- 
mosphere while some are trapped in the corona 
after being accelerated. The physics of how ener- 
getic electrons propagate from the corona to the 
chromosphere, and the contrary, how confinement 
of electrons happens, remains poorly understood. 
The transport of electron heat flux has previ- 
ously been modeled as conduction mediated by 



classical collisions (Spitzer 1962), convection at a 



subthermal characteristic speed observed in simu- 



lations of laser fusion studies ( Manheimer & Klein 



1975), and turbulent transport limited by anoma- 



lous resistivity ( |Manheimer|1977[|Tsytovich|1971j ), 
also called anomalous conduction. Recent analyti- 
cal and numerical studies of the transport of super- 
hot electron fluxes with energies greater than 10 
keV show that the classical conduction model is in- 
applicable because the mean-free-path of electrons 
exceeds the characteristic length scales of the sys- 
tem ( Oreshina fc Somov||2011 ). 

In the turbulent transport scenario, anomalous 
resistivity arises from electron scattering by tur- 
bulent wave fields excited by instabilities involving 
the interaction of the energetic electrons, ambient 
electrons and/or ions. One commonly considered 



mode is the ion-acoustic instability (Manheimer 
T977i |Tsytovich|[T97T| [Smith fc Limequist||1979| ) 



When accelerated electrons propagate outwards, 
ambient electrons are drawn in as a return cur- 
rent. Their relative drift with ions excites ion- 
acoustic waves, which scatter the energetic elec- 



trons and anomalously enhance resistivity (Man- 



heimer 1977). The heat flux carried by the ener- 



getic electrons is limited by the saturation of the 
instability due to ion trapping. 

The formation of an anomalous conduction 
front due to ion acoustic turbulence was consid- 
ered as a means to confine hot electrons for the 



production of hard X-rays (Smith & Lilliequist 
1979). The model, based on a ID one- fluid code 



with a grid spacing of 100 km (~10'* ion inertial 
lengths), could not capture processes occurring at 
electron scales. However, it was argued that the 
front propagates at approximately the ion acoustic 
speed Cg and, with some heuristic analysis, pre- 
dicted that the front would become thinner than 
the grid spacing, although this was artificially 
disallowed to avoid numerical instability. More 



recently, the existence of a thermal front, literally 
defined as a region that links plasmas in thermal 
nonequilibrium and sustains the temperature dif- 
ference for longer than the electron free-streaming 
time, was studied in ID electrostatic Vlasov sim- 
ulations ( [Arber fc Melnikov|2009[ ). The main goal 
was to refute the results of an earlier particle- 



in-cell (PIC) simulation (McKean et al. 19901 in 



which a conduction front was not seen. The au- 
thors observed the formation of a temperature dif- 
ference propagating at a speed comparable to Cs 
and identified the behavior as a conduction front. 
The physics of the responsible mechanism was not 
identified or investigated. We suspect that front 
to be an ion acoustic shock and the temperature 
difference to be a result of shock heating due to 
its extremely sharp transition. We observe simi- 
lar behavior in our simulation, to be presented in 
section |3l 

To gain a better understanding of the problem 
of energetic electron transport from the corona 
to the chromosphere, an important question that 
must be answered is whether thermal conduc- 
tion is suppressed when the accelerated electrons 
start to propagate into the immediate surrounding 
plasma. Our goal in this paper is to answer this 
basic question. 

To study the physics at the electron scale and 
the nonlinear aspects of the problem, we perform a 
2D PIC simulation. We observe that the transport 
of energetic electrons is suppressed by a double 
layer that occurs at the shortest scale in a plasma, 
the Debye length Xue- 

A double layer (DL) is a nonlinear electrostatic 
electric field localized within two adjacent layers of 



equal and opposite net charge (Block|1978 Raadu 



fc Rasmussen|1988[ [Singh et al.|1987i . It is highly 
localized to ~10 Xoe because quasi-neutrality in 
a plasma can only be violated on Xoe scales. It 
can be understood as a capacitor in a plasma. An 
ideal DL is a unipolar electric field, which corre- 
sponds to a monotonic drop in the potential across 
the structure. In most situations, however, the po- 
tential can contain dips and bumps at the high or 
low potential sides, but with a net potential drop 
across the entire layer. The potential drop (j) can 
accelerate, decelerate and reflect particles enter- 
ing from the two sides of the DL. For example, 
electrons incoming from the high side will be 
decelerated (or reflected) if their kinetic energy is 



greater (or less) than (j). DLs can be classified into 
strong and weak depending on whether (p is much 
greater or comparable to the mean energy of the 



reflected particles on either side of the DL ( Raadu 



& Rasmussen 19881. Much research interest on 



DLs has focused on their ability to directly accel- 
erate particles in auroral (ionospheric and magne- 
tospheric) and astrophysical plasmas ( Singh et al. 



The simulation domain is of course far smaller 
than any realistic flare loop. 

The electrons are initially modeled by bi- 
Maxwellian distributions and the ions by a Maxwellian 
distribution. At the simulation's start, one-third 
of the electrons, centered at the middle of the 
domain, have temperatures parallel to the back- 



ground magnetic field T), n/r 



''I'-A' 



=0.5 (orange in 



1987|). For particle acceleration, the potential Fig. [T]) while the remaining two-thirds have 



across the DL is usually maintained by some ex- 
ternal energy source. Mechanisms that have been 
used in simulations to produce DLs include an im- 



posed potential drop or current ( Singh et al. 1987 1 . 
In the latter case, DLs occur as nonlinear waves 
arising from current-driven instabilities. The DL 
observed in our simulation is generated by impos- 
ing a large field-aligned temperature jump in the 
initial state which results in strong currents that 
drive the DL. 

Parameters suitable for flare settings are used 
in our simulations and described in section [21 In 
section [Sj we present the results. A brief discus- 
sion of our simulations and application to con- 
fining above-the-looptop hard X-ray sources are 
given in section HI We summarize in section [5] 

2. SIMULATION 



We are exploring the problem of energetic 
electron transport from the corona to the chro- 
mosphere by a two-dimensional electromagnetic 
particle-in-cell simulation using the p3d code 



(Zeiler et al. 2002). The energetic electrons are 



pre- heated (i.e., present at the start of the simu- 
lation). We consider it to be more natural that 
the hot population not have a preferred direction 
of propagation, so the simulated energetic elec- 
trons are not beamed in the initial state. Fig. [T] 
is a cartoon of the initial simulation setup. Our 
system represents a symmetric local segment of 
a fiare loop with very hot electrons centered at 
the loop-top in contact with ambient cold elec- 
trons. The computational size is L^ x Ly= 655.36 
X 2.56 dl with a cell size of 0.02x0.02 d^. There 
are 400 particles per cell, x is the direction par- 
allel to the initial background magnetic field Bq. 
de=c/uipe is the electron inertial length, where 
Wpe=(47rnoe^/rne)^'^ is the electron plasma fre- 
quency and no the initial density, which is uniform. 
The boundaries are periodic in both dimensions. 



Tg ||/TO.iC^=0.1, corresponding to ratios of parallel 
plasma pressure to magnetic pressure of /3e||=l 
and 0.2, respectively. CA=i3o/(47rTOino)^/^ is the 
Alfven speed based on Bq. The temperature per- 
pendicular to Bq is Tf,^_\_=Q.l=Ti everywhere, so 
/3e__L=0.2=/3i. Since it is easier to accelerate elec- 
trons along the field lines, the hot electrons are 
taken to be anisotropic. Ambient electrons and 
ions are taken to be isotropic. Note that the unity 
/? hot electron pressure is consistent with recent 



over-the-limb fiare observations (Krucker et al 



2010) 



Other parameters used are: mass ratio 77ie/TOi=l/100 
and speed of light c/cyi=100. Electric fields are 
normalized to Eq—caBq/c. Thermal speeds are 
defined as Vix — C^Tx/mxy^^ . The Debye length 
based on the initial hot electron temperature is 

^De—VteO,h/^pe—0-^ de- Wc will USC de aS the 

unit of length, but it can be conveniently con- 
verted to Xdc- We will use T,, to represent Tg y in 
the following for simplicity, since it is the parallel 
heat transport that is of most interest. The total 
duration of the simulation is less than the transit 
time of the domain by electrons at 1.5 Vteo,h, so 
the majority of hot electrons will not encounter a 
boundary during a run. 

3. RESULTS 

3.1. A jump in electron temperature 

The primary drive of the evolution of the sys- 
tem is the contact of hot electrons with cold elec- 



Ly=2.56 d, 

_ + 



Fig. 1. — A cartoon of the initial simulation setup. 



trons, so the dynamics starts in the region of con- 
tact. Since this is a symmetric system, it suffices 
to study either side of the contact. We choose 
the right side, where hot electrons propagate out- 
wards with positive velocities. The figures are av- 
eraged in y since there is no significant variation 
in that direction. Fig. [2] shows the time evolu- 
tion of electron temperature in the region where 
the hot and cold electrons come into contact. The 
initial profile (black) represents a high tempera- 
ture region decreasing to a low temperature region 
through a narrow transition that has a length scale 
of several de's. We do not expect the results to 
be sensitive to this transition scale since we find 
that the DL develops from a streaming instabil- 
ity after the temperature discontinuity has broad- 
ened. Over a short period of time ('^lOO u}~^ , 
blue) , the transition broadens as the hot electrons 
free-stream into the cold electron region. If free- 
streaming were to continue, we would expect the 
profile to continue to broaden due to the mixing 
of the two populations. However, a small jump 
develops (see a;pet=300 ) within the smooth tran- 
sition, meaning that the escape of hot electrons 
by free-streaming is inhibited. Over time, the 
jump grows in size and significantly suppresses 
mixing due to free-streaming. As a result, the 
previously smooth transition is converted into a 
distinctive jump across two separate temperature 
regions. This indicates that the transport of hot 
electrons into the surrounding plasma is being sup- 
pressed. Other features of Fig.[2]include Langmuir 
waves propagating in the cold electron side (see 
e.g., the wiggles on the right side of the curve at 
Wpet^SOO) that are excited by the bump-on-tail 
instability as hot electrons free-stream into that 
region. Beam modes are excited on the hot elec- 
tron side (see, e.g., longer scale waves on top of 
the curve at Wpet=1400) by the electron-electron 
streaming instability as cold electron beams that 
are accelerated by a DL enter the hot side and 
interact with the hot electrons. 




Fig. 2. — Evolution of electron temperature Te 
parallel to the local magnetic field near the tran- 
sition between hot and cold electrons, at x= Odg. 
We begin with two electrons populations in con- 
tact (black) and let them evolve. The initial gradi- 
ent becomes a smooth diffusive profile at Wj,et=100 
(blue). The smooth transition is then interrupted 
by a jump which signifies transport suppression 
from the high to low Tg sides. 
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3.2. DL evolution and strength 

The heat fiux suppression comes from a DL that 



arises within the transition region. Figure 3(a) 
shows the time history of the electric field parallel 
to the local magnetic field, E^, and |3(b)] the time 
evolution of the electric potential jump across the 
DL normalized to the instantaneous hot electron 



Fig. 4. — e(j)/Te^h (black) and Tg (green), overlaid 
with E^ (blue) and the zero E^ position (dotted 
cyan), at a;pet=1200. The DL is characterized by 
a monotonic drop in the potential. The potential 
barrier suppresses the transport of hot electrons 
and results in a drop in Tg across the DL. 
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Fig. 3. — Evolution of (a) E^, the parallel electric field and (b) e(j)DL/Te,h (solid black), the electric potential 
jump across the DL normalized to the instantaneous hot electron temperature at the center of the hot region 
(dotted line is from a second simulation with Teo,/t=l that will be discussed later), overlaid with Nref/No 
(solid green), the fraction of return current electrons reflected at the foot of the DL divided by the total 
initial electron number. 



temperature measured at the center of the hot 
region, e(t>DLlTe,h- Te^h in e(t>DL/Te^h thus rep- 
resents the core temperature of the hot electron 
population. T^.h drops over time as can be seen 
on the left side of Fig. [2] because hot electrons 
continuously leak out. We normalize the poten- 
tial jump of the DL to T^.u to more accurately 
reflect the suppression strength of the DL on the 
hot electrons that remain in the source region at 
any particular time. Also shown in Fig. |3(b)| is 
the fraction of return current electrons that are 
reflected at the foot of the DL, normalized to the 
initial total electron number, Nref/No. The reflec- 
tion mechanism and the reason for its importance 
are discussed later. E^ is averaged over an electron 
plasma period to eliminate initial fluctuations in 
the contact region that decay over time. The DL 
emerges at Wpet~100 around x=0. At Wpet~550, 
it emits a large amplitude shock that moves in the 
positive X direction. The DL is the more struc- 



tured dark region in Fig. |3(a) that slowly drifts 



in the negative x direction. The DL is constantly 
evolving over the entire course of the simulation. 
At a;pet=1200, it reaches its maximum amplitude 
and extends over a scale of 4 dg, centered at x~ 
-6de- Fig. |4]is a cut of Fig. 3(a) at this time in a 
blowup around its location. Shown are E^ (blue), 
the electric potential normalized to the core hot 
electron temperature e4>/Te^h (black) and the elec- 
tron temperature Tg (green). E^ has a large pos- 
itive peak at x~ -Qde that causes a large poten- 
tial drop that reflects hot electrons and produces 
a sharp drop in the electron temperature. We de- 
fine the strength of the DL by its potential jump 
e0DL which is calculated as follows. To the left 
of the DL in Fig. |4j the potential increases, oscil- 
lates and eventually reaches a fairly constant value 
(near x~ -13de), the high end value. To the right 
of the DL, the potential drops and oscillates about 
a roughly constant value (near x~ -5de), the low 
end value. e4>r)L is the difference between the high 
and low end values and is shown as a function of 
time (solid black) in Fig. 



3(b) e4>DL/Te,h 



measure of the suppression capability of the DL 
on the hot electron transport. If the DL is large 
enough such that e<j)DL/Te^h ~ Ij the thermal bulk 
of the hot electrons will be reflected by the poten- 
tial barrier of the DL. This will imply a significant 
reduction of the electron heat fiux as well. 

Note also in Fig. |4]that the DL is not merely a 



unipolar electric field. It has a small negative leg 
between x= -6 and -Sdg, on the low Tg side. While 
the bigger positive part suppresses streaming hot 
electrons coming from the left, the negative leg can 
reflect cold electrons coming from the right (low 
Te side). The reflection of these electrons is the 
drive mechanism of the DL, which we now discuss. 
Fig. Is] is a cut of Fig. 3(a) at CL;pet=150, shortly 
after the DL emergence. Multiple peaks in E^ de- 
velop in the Tg transition region. In the linear 
regime (a;pet<150), these waves have small ampli- 
tude, comparable to the background fluctuations. 
They gradually grow into large-amplitude nonlin- 
ear structures at later times. The peak which has 
the largest amplitude dominates the others and 
becomes the DL. It then starts to hinder hot elec- 
tron free-streaming and produces an obvious drop 
in Te as can be seen at a;pet=150. The two smaller 
peaks on either side of x=0 later turn into fairly 
symmetric waves, which produce no significant po- 
tential jump. The wave turbulence from which the 
DL emerges is driven by the return current inter- 
acting with the background ions. The DL is also 
driven by this streaming instability, whose origin 
is now described. 



3.3. 



DL drive mechanism: streaming in- 
stability 



Figure 6(a) shows the electron phase space at 



a;pet=150. The hot and cold electron distributions 
are at the left and right sides, respectively. Note 
the narrower thermal spread of the cold electrons 
at the right, as well as the wider spread of hot elec- 
trons at the left. The initial temperature transi- 
tion is near x=0. A beam of cold electrons is enter- 
ing the hot region. At uipet—Q, both distributions 
are non-drifting and there is no beam. As the hot 
electrons free-stream into the cold electron region, 
ions cannot follow and are left behind as a charge 
imbalance. The resulting electric field drags a cold 
electron beam towards the hot electron region as 
a return current. The relative drift between this 
return current and the ions excites the instability. 
Fig. |6(b)| is a cut of the phase space at the po- 
sition of the DL (x~ -Ide). The return current, 
with a negative beam speed, sits on top of the 
streaming hot distribution, which is moving to the 
right and has a wider thermal spread. The drift 
of the return current electrons with respect to the 
ions (red) excites the ion/return-current-electron 
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Fig. 6. — (a) Electron phase space at ujpet=150 and (b) a cut of (a) at the position of the DL at x^ 
-Ide, indicated by a vertical black line. The ion distribution fi (red) is overlaid on the electron distribution 
fe (black). Since the ions and cold electrons have the same temperatures, the ion thermal speed is ~ 
■\/mi/nii,=10 times lower than the return current thermal speed. The integrated areas under the curves 
(i.e., their densities) are the same, so fi peaks at ^10 times the height of /g. The vertical scale is normalized 
to the maximum of fi. 
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Fig. 5. — The electric field E^ at ajpet=150 , over- 
laid with Te (green). The large amplitude wave 
that forms around x<0 arises from an ion/return- 
current-electron streaming instability. 



streaming instability. The linear growth rate 7 of 
the ion-electron streaming instability (also known 
as the Buneman instability (IBuneman 1958)) for 



the most unstable wave calculated in fluid theory 
is given by 



^/3 / me 
^ \2mi 



1/3 



(1) 



Using the parameters of our simulation, 7 ^ 0.15 



'.'pe- 



The growth time r is thus r ^ 7 



-7 Wp/. 



The instability is sufficiently strong to produce the 
DL at Wpet-^IOO. 

The time evolution of the electron phase space 
reveals that the growth of the DL depends on the 
number of return current (RC) electrons reflected 
at the foot of the DL. One can see from Fig. |4] 
that the electric field Ex of the DL has a nega- 
tive leg on the low Tg side. This shows up as a 
potential dip on the low potential side. The RC 
electrons are reflected by the negative leg, which 
effectively reduces the RC speed and thus con- 
verts the kinetic energy of the RC to the electro- 
static energy of the DL. By exchanging momentum 
with the reflected electrons, the DL gets amplified 
(Lotko 1983). Electron reflection also indicates 
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Fig. 7. — Electron phase spaces at: (a) ujpet—550; and (b) ujpet=700. The vertical black lines mark the low 
potential side of the DL and to its right is the negative leg E^ of the DL that can reflect the return current 
electron beam. 



when the instability enters into a nonlinear regime 
(a;pet>150). The DL growth depends on the num- 
ber of RC electrons reflected at the foot of the DL. 
Qualitatively, the dependence of DL growth on the 
number of reflected electrons can be seen by con- 
trasting the phase space at a time when the DL 
is growing, e.g., at Wpet=550, with a time when 
it is not growing, e.g., at Wpet=700 (as e(j)DL/Te.h 
plateaus in Fig. |3(b)| ). At Wpet=700 (Fig. |7(b)[ ) 
to the right of the DL at x> -4de (vertical black 
line), there are many fewer reflected electrons than 
to the right of the DL at a;pet=550 (Fig. 7(a) ), x> 
-2de. We estimated the number of reflected RC 
electrons, Nref, by integrating the electron distri- 
bution to the right of the DL over a distance of 
one de, which is the typical width of the DL, and 
for Vex G(0, 0.4 Vte,h), which is approximately the 
velocity range covered by the reflected RC elec- 
trons as they turn around and acquire positive 
velocities. N^ef is normalized to A'o, the initial 
electron number contained in the same interval. 
Note that the integration for Nref also includes 
a contribution from hot free-streaming population 
in that velocity range. However, this contribution 
stays nearly constant and the variation in N^ef 
mainly comes from the change in the number of re- 
flected cold RC electrons. Therefore, Nref suffices 
as an estimate of the reflected RC electrons (plus 
some constant offset). Nref/No (green) is plot- 



ted on top of the DL size, ecj^DL/Teji, in Fig. 3(b) 
to show their correlation. As Nref/No increases 
during Wpet<550 and >950, ecffoL/Te^h bas a gen- 
eral increasing tendency as well. As Nref drops 
abruptly at Wpet=550 and stays at a low level un- 
til Wpet=950, e(j)DL/Te,h plateaus. Hence, there is 
good correlation between the DL growth and the 
number of the reflected RC electrons. 

An examination of the evolution of the electric 



field in Fig. 3(a) reveals that at Wpet^SSO, the DL 
sheds an ion-acoustic shock wave, which shows up 
as an ion hole in ion phase space. Once it leaves 
the DL, the shock accelerates to the right while the 
DL continues to slowly move to the left. There- 
fore, their separation increases. The shock wave 
is bipolar and reflects most of the RC electrons 
with its large-amplitude negative leg. One can see 
in Fig. 7(b) at Wpet=700, most of the RC elec- 
trons are reflected behind the shock at x< -2de- 
Nref /No shows an abrupt drop at ujpet—?)^Q be- 
cause most of the RC electrons are stuck behind 
the shock and do not reach the DL. The shock dis- 
sipates through ion heating. The ion temperature 
increases by as much as a factor of two along the 
path traversed by the shock. The shock gradually 
diminishes in amplitude. Similar dynamics (accel- 
eration and decay) of an ion hole mode at the low 
potential side of a weak DL created in ion-acoustic 
turbulence is also observed in current-driven sys- 



terns ( Barnes et al.||1985 ). 

During our simulation, e(f>DL/Te.h reaches 0.8 
at Wpet=1200 (Fig. |4|. This means that hot elec- 



trons with velocities v/vte,h < v 0.8 ~ 0.9 will be 
reflected. The electron phase space at ujpet—1200 
in Fig. [8] reveals that the population of hot elec- 
trons with v/vteji < 0.9 wraps around from pos- 
itive to negative velocities when encountering the 
DL. So they are reflected and hence confined in 
the source region. The suppression observed here 
is significant. We calculate the heat flux carried 
by kinetic hot electrons that is po sitively directed 

oo 

as Q{x,t)—^mf. J v^fe{x,v,t)dv (Krall & Trivel- 





piece 1973). The mean flow speed is <2% of the 



hot electron thermal speed and hence negligible in 
the calculation of Q{x,t). Similar to Tg, Q has a 
sharp jump across the DL. Towards the end of the 
simulation, Q is suppressed by ~40% to the right 
of the DL compared with the higher value to the 
left of the DL. 

3.4. Dependence of DL strength on hot 
electron temperature 

To determine the dependence of the DL po- 
tential drop on the hot electron temperature, we 
performed another run in which we raised the 
initial hot electron temperature T'e,/i.|t from 0.5 
to 1, and Bq from 1 to a/2 such that /3e.h,|| re- 
mains unity, while keeping all other parameters 
unchanged. Slightly increasing Bq should not have 
an effect on the result because the phenomena we 
observed here are dominantly electrostatic. 

The size of the DL depends on the hot electron 
temperature. The DL potential drop e(t)DL/Teji 
from another simulation with a higher initial hot 
electron temperature T'eO,/i=l is plotted as a dot- 
ted line in Fig. |3(b)| This potential drop is nor- 
malized to the instantaneous hot electron temper- 
ature at the center of this new hot region. The DL 
is on average stronger than that in the simulation 
with initial Teo,h=0.5 (solid line). The basic dy- 
namics of the DLs are the same in both cases. A 
phase of growth is followed by a period of pause, 
and the whole cycle repeats. The DL evolves and 
grows faster in the higher T^o^h run, probably be- 
cause the return current electron response, which 
is the driver of the DL, is faster due to the in- 
creased thermal speed of the hot electrons (relative 



to those of the ambient plasma) in the Teo.h=l run 
than in the Teo.?i=0.5 run. There is also a mod- 
est increase in the suppression of heat flux Q{x,t) 
from -40% in the Teo./i^O.S run to ~50% in the 
higher Teo^h run. Therefore, one might expect a 
higher ratio of hot to cold electron temperatures 
would lead to a stronger DL and associated sup- 
pression of transport. 

4. DISCUSSION AND APPLICATION 

We have studied the suppression of electron 
thermal transport due to double layer (DL) for- 
mation in 2D PIC simulations. However, our sys- 
tem is effectively ID since the transverse direction 
is only « 20A£)e wide. The presence of the sec- 
ond dimension enables us to average the data for 
noise reduction. 2D effects on DLs are not investi- 
gated in this study, but do not seem to be signifi- 
cant. Extended (planar) weak DLs over 80 Xd,, in 
transverse extent were observed in earlier 2D PIC 
simulations where subthermal electron current was 
injected at the boundaries of the parallel direction 
(Barnes et al. 1985). Such electron injection may 



correspond to the return current electrons being 
drawn in at the contact region of the two electron 
populations in the present simulations. The pla- 
nar 2D DLs tend to develop substantial substruc- 
ture across the magnetic field late in time, but the 
basic dynamics of the DLs are the same as in ID. 
We expect the formation of DLs, and therefore 
transport suppression of hot electrons, to persist 
in 2D. However, 2D DLs are unlikely to be planar 
over ~ lO^Aoe, the scale of coronal looptops. On 
the other hand, a variation of the DL over a large 
transverse distance is not likely to strongly impact 
the largely ID dynamics studied here. The injec- 
tion of a narrow band of hot electrons, however, 
might behave very differently. This case remains 
to be studied. 

We also note that the computational size of 
the PIC simulations, which is 655(ie — 66m in 
the corona, is small compared to the actual size 
of flaring loops. They, however, reveal important 
physics down to the kinetic electron scale, the De- 
bye length, where DLs form. We expect the DL 
to continue strengthening in a larger domain that 
allows the system to evolve for a longer period. 

Applying our results to the above-the-looptop 



hard X-ray sources (Krucker et al. 2010 20071 that 



show an exponential decay on a timescale of more 
than two orders of magnitude longer than the elec- 
tron transit time through the source, the forma- 
tion of a DL provides a plausible mechanism to 
confine energetic electrons in the source region. 
Electrons that have kinetic energies less than the 
potential drop of the DL are reflected back to the 
source region. Over the entire course of the simu- 
lations, the DL tends to amplify through multiple 
phases of growth. There are no signs of decay. 
Since it is driven by the return current electrons, 
which exist for as long as some hot electrons are 
leaking out, the DL can persist until the eventual 
depletion of all hot electrons in the source region, 
which should take much longer than the electron 
transit time through the source due to the strong 
transport suppression caused by the DL. The DL 
lifetime is therefore expected to be much longer 
than the transit time. This can naturally explain 
the prolonged lifetimes (>100 s) of the energetic 
electrons from the above-the-looptop hard X-ray 
observations compared with the electron transit 
time (~1 s) through the source. 

5. CONCLUSION 

Results from PIC simulations of a pre-accelerated 
hot electron source in contact with ambient 
plasma, which provide a basic model of coronal 
looptop sources, are presented. The hot electrons 
free-stream along the magnetic field initially, but 
transport is suppressed by the formation of a dou- 
ble layer that is driven by an ion/return-current- 
electron streaming instability due to the drift of 
return current electrons with respect to the ions. 
The DL subsequently reflects and decelerates hot 
electrons as they begin to propagate into the sur- 
rounding plasma. A significant reduction in elec- 
tron heat fiux through the ambient electron region 
is observed. Such a suppression mechanism differs 
from the conventional picture of electron scatter- 
ing by turbulence. Our results are consistent with 
observations of both free-streaming propagation 



(Aschwanden et al. 1996) and the suggested con 



finement at looptop (Masuda et al. 1994 Krucker 



et al. 2010 2007) of hard X-ray producing elec- 



trons, corresponding to the escaping and reflected 
hot electrons, respectively. 
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Fig. 8. — Electron phase space at Wpet=1200 with 
e4>DL/Te,h ~0.8. To the left of the DL location, 
x~ -6de, hot electron populations with v/vte,h < 
0.9 (cyan and green in color) are reflected by the 
potential barrier. They are trapped in the hot 
region. For higher velocity electrons (at the tail; 
blue in color), there is a drop in velocity as they 
are decelerated. Both contribute to suppression of 
the heat flux. 
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